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Abstract. Magnetization dynamics in nanomagnets has attracted broad interest 
since it was predicted that a dc-current flowing through a thin magnetic layer can 
create spin-wave excitations. These excitations are due to spin-momentum transfer, a 
transfer of spin angular momentum between conduction electrons and the background 
magnetization, that enables new types of information processing. Here we show 
how arrays of spin-torque nano-oscillators (STNO) can create propagating spin-wave 
interference patterns of use for memory and computation. Memristic transponders 
distributed on the thin film respond to threshold tunnel magnetoresistance (TMR) 
values thereby detecting the spin-waves and creating new excitation patterns. We show 
how groups of transponders create resonant (reverberating) spin-wave interference 
patterns that may be used for polychronous wave computation of arithmetic and 
boolean functions and information storage. 
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1. Introduction 

Polychronous Wavefront Computation (PWC) as a paradigm for a post digital era 
was introduced by Izhikevich and Hoppensteadt in 2009 pQ following the early work 
of Izhikevich [2], where the importance of spatial propagation and axonal delays 
in the brain were investigated. The key ingredients for PWC are a medium that 
can support interference patterns of propagating activity packets, and transponders 
that can sense the size of incident activity and respond to super-threshold inputs by 
generating a propagating wave. Mathematical models that support such features are 
two dimensional wave equations and Schrodinger's equation in bounded domains with 
absorbing boundary conditions. Transponders are devices located at specific locations 
that can receive and transmit activity when appropriately stimulated. PWC is inspired 
by studies of the brain by Izhikevich [2j, where it was shown that delays caused by 
axonal conduction times introduce new possibilities for network behavior, including 
memory and computation. These are referred to as being polychronous aspects of brain 
function. 

Here we describe how to construct devices that use polychronization arrays of 
spin-torque nano-oscillators (STNOs) in two-dimension ferromagnetic thin films. The 
propagation times of peak wave interference correspond to conduction-time delays. The 
analogy of excitation of a neuron is excitation of a transponder that may fire when 
a certain threshold of inputs from other transponders is exceeded. In this case, the 
transponder generates a new wave, perhaps after some time delay. The work presented 
here is not an attempt to model a brain, but to accomplish some brain-like behaviors 
on much shorter time scales and on smaller spatial scales. Our devices support the 
propagation of activity patterns on spatial scales of nanometers at gigahertz frequencies, 
while brain activity occurs on millimeter scales at tens-of-hertz frequencies. 

1.1. Wavefront Intersections 

Associated with a wave is a cone in space-time whose cross-section expands in time 
proportional to the wave speed. If the region is isotropic, the cones are right cones. If the 
region is irregular, the cones are irregular, as well. Energy packets may not be uniformly 
distributed on the cone, as we demonstrate later for spin waves. The surface of a cone 
will be referred to as being a wave front along which energy packets may propagate. 
Interference patterns created by energy packets on cones reveal the possibilities for 
PWC and guides programming. In some applications, energy is distributed uniformly 
on the cone (e.g., the excitation produced by a single STNO diffuses isotropically) . 
In others, energy packets, while always moving on the surface of the cone, may not 
be uniformly distributed on it (e.g., certain excitations of some STNO arrangements 
produce interference patterns and consequently irregular energy distributions by means 
the space-time cones). 

The notion of signal conduction delay has a clear representation in interference 
patterns created by waves. Two transponders firing at different times may excite 
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propagating energy packets that may be visualized as lying on cones in space-time. 
A cone will begin at the location and time of initial excitation, and expand outward 
with increasing time. Two cones may intersect in a conic section, usually in a parabolic 
curve. This curve describes the evolution of the intersection of the two waves. The 
parabolic-shaped curve on which wavefronts intersect will be referred to here as being 
an intersection parabola as shown in Fig. [TJ (If one cone is contained within the other, 
there will be no intersection.) A transponder lying on an intersection parabola may be 
excited by the wavefronts converging on it. Transponders may be placed at differing 
distances from a common target transponder. Given that each wavefront has the same 
finite propagation times, delays are encoded in Euclidean distances. We also allow 
transponders to be excited externally. We can therefore introduce additional time delays 
into our system simply by exciting two transponders at different times. 

Complicated polychronous activity is coded in PWC by the relative positions of 
transponders, and in the timing of external excitations. We illustrate this idea in Fig. 
Q] by describing possible curves of intersection of pulses emanating from two separate 
transponders Ti and T 2 generated at separate times, say r units apart. 
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Figure 1. Transponders T\ and T2 fire with a range of time delays. The waves will 
intersect along the curves indicated in dashed lines projected in to the t = plane. 
Transponders labeled with stars detect and eventually fire whenever there is an inter- 
pulse interval of length Ti, respectively; the array processes signals from transponders 
Ti and Ti by converting temporal delays into spatial firing patterns. 

Several aspects of programming for PWC are comparable to sound location and 
response mechanisms in animal brain. For a given task, we show how to arrange 
transponders to accomplish it. This mode of programming is comparable to learning 
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in the human brain. For example, multiplication tables are learned by shaping our 
neural network to provide a look-up table mechanisms for the result of multiplying two 
numbers. One aspect of the work here describes how to create look-up table arrays 
using PWC. Other aspects regard storage of computations, detection of inter-wavefront 
intervals as a means for recovering stored information, and possible physical mechanisms 
for implementation. 

It is shown in [H [3] how to place transponders to accomplish computations such as 
addition, subtraction, multiplication and division. 



2. STNO design 



The interaction of a spin-polarized current with a magnetic film results in magneto- 
resistance |ll |5] and a spin-transfer-torque effect [6l [7J El [9j [10]. The latter consists 
of a torque exerted by the spin-polarized electrons on the background magnetization. 
Promising applications based on the spin transfer effect are spin-torque nanoscillators 
(STNOs) [HI QUI E] where spin-transfer-torque is used to drive a GHz oscillation of the 
magnetization direction of the magnetic film. STNO consisting of a point contact to a 
thin film ferromagnet (FM), were first proposed theoretically in 1996 [El [7]. DC current 
densities (j > jc) generate a high-frequency dynamic response (1 to 100 GHz) in the 
FM layer and result in the emission of spin- waves [8]. 
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Figure 2. Schematic of an STNO. (a) A single STNO: the charge flows through a 
point contact to a thin ferromagnetic layer. The magnetic moments in the contact area 
precess, exciting spin-waves that diffuse in the film, (b) A double STNO arrangement 
is shown. The double point contact creates a spin-wave interference pattern, (c) The 
spatial dependence of the amplitude of the excitation |m| 2 = m 2 + rriy is plotted as 
well as one of the in- plane components, m x . (d) Three point-contacts in a triangular 
arrangement are excited with a current pulse producing a magnetic excitation; the 
magnetic excitations diffuse and interfere throughout the two-dimensional thin film. 
Intersecting waves create six wave-packets that propagate in different directions. 
Simulations are performed in a 500 by 500 nm film with point-contacts of 20 nm 
diameter, and separation between contacts of 60 nm. 
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The basic concept of STNO devices is illustrated in Fig. [2(a). A nanoscale electrical 
contact is attached to a multilayered ferromagnetic structure [9J QI)J [TT] . The multilayer 
consists of a fixed and a free magnetic layer separated by a non-magnetic metal. The 
fixed layer determines the polarization of the current that excites magnetic moments in 
the free layer. A typical fixed layer would be a 10 nm thick permalloy film (with in-plane 
magnetization) or a multilayered structure (with an out-of-plane magnetization). The 
free layer may be a much thinner film (1-3 nm) [T21 The contact sizes range from 
tens to hundreds of nanometers. Phase locking between oscillators has been observed 
for distances up to 500 nanometers [UJ [151 [IS]> an d the spin waves are predicted to 
propagate distances of tens of microns. 

The magnetic moments of the free layer initially precess around the direction of the 
applied field and eventually align with the field as the result of damping. The response to 
a (polarized) current in a free magnetic thin layer is phenomenologically described with 
an additional current- dependent positive damping term in the Landau-Lifshitz-Gilbert 
equation. The result is a microwave excitation of magnetic moments of the thin film 
layer. Such excitations have a precession frequency due to the local magnetic field that 
depends non-linearly on the amplitude of the excitation. If we consider a unit vector to 
represent the direction of magnetization, rh, then in the geometry considered in Fig. [2] 
the out-of-plane magnetic moment m z does not oscillate and it represents the envelope 
of the modulated in-plane components, m x and m y , (m 2 z = 1 — (m x + m 2 ) = 1 — m 2 , 
we call m = (m x ,m y ) the in-plane component for convenience). Weak excitations at a 
single contact decay radially (see, Fig. EJ^c)). 

Micromagnetic modeling on shaped magnetic thin-film waveguides [17] and on 
double coupled STNOs [18] has been done, and interference patterns have been found 
[T9| [20] . Arrangements of more than one STNO can create interference and enhance 
activity in some preferred directions in the film's plane. The richness of the magnetic- 
moment dynamics during the dissipation of initial pulses is indeed enormous and 
it motivated our methodology for generating interference patterns throughout the 
magnetic layer. A spatial pattern created by three STNOs arranged in the vertices 
of a triangle is depicted in Fig. [2(d). 

Multiple current pulses excite multiple spin-waves that interfere. These create 
patterns in the envelope of the magnetic moments (i.e., m z for our geometry), and the 
constructive interference of wave-packets from several STNOs creates peaks in particular 
radial directions. By engineering both the distance between contacts and their sizes, it 
is possible to control spin-wave patterns. Notice that excitations must be phase locked 
within the STNO arrangements; an external background signal might be used to lock 

in phase or in frequency the STNOs m[igi2D|221[2ai21[251l2Ell2Zl. In addition ' 
the frequency modes in the background signal might contribute energy for information 
storage [28] (see our discussion in section Appendix A). 
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3. Mathematical Theory 

Mathematical modeling is based on the Landau-Lifshitz-Gilbert equation in thin 
ferromagnetic films with an additional positive damping term corresponding to the 
transferred torque from the polarized current to the film. Let M denote the 
magnetization vector in the thin film, then if an external magnetic field is applied, 
the spin-magnetic moments will initially precess around the direction of the applied 
field, and eventually align with the field as the result of damping. An applied polarized 
current will cause a positive damping, and if the applied current exceeds a certain 
threshold (j > jc), then a spin- wave excitation will be created [6j [TJ El El [TO]. This is 
described by the equation 
dM 

— = -ItI^oM x H eff — ctM x (M x H eff ) + /3(x)(M x M x m p ), (1) 

where the precession (first term) and damping (second term) are controlled by the 
effective field H e ff , being the sum of the applied field, the demagnetizing field, and the 
exchange field, 

H eff = H z - M z z + —JL—V 2 M. (2) 

\-f\/j,oM s h 

The spin-torque (third term) is controlled by the spin polarization direction of the 
applied current, m p . No variations in z are considered (the Laplacian is taken in two 
dimensions), as the free layer is thin with respect to the exchange length. The function 
/3(x) is a Heaviside function defining the contact sizes and locations. Furthermore, the 
function /3 depends on the current intensity, the layer thickness and the spin polarization. 

We define dimensionless parameters 

M uj m , x 
m = — — = (m x ,m v ,m z ), t = r, x = - — , (6) 

where um = ^1^qM s is the Larmor frequency for an applied field M s and Z e x — 



D/"ffj,oM s h is the exchange length. 

Since the magnetization vector m lies on the unit sphere and the motions of interest 
are small deviations from equilibrium, we consider the components of m as 



m = m x + vm y , m z — y 1 — |m| , 

With these variables the Landau-Lifshitz equation takes the form of a nonlinear 
Schrodinger equation 
dm 

i—— = —mV 2 m z + (m z + ia') V 2 m — (m z — h)m — ictih — m z )m z m + ia' (I VmJ 2 + I Vml 2 ) m 
or 

(m z m)m Fz + (m x m y - (m 2 z + ■m 2 y ))m Fx + (m x m y - (m 2 z + m 2 x ))m Fy 

din 

— Im(m*V 2 m) + a'(h — m z )\m\ 2 + a' [V 2 m z + (|Vm| 2 + \Vm z \ 2 )m z ] 



(4) 



Or 



\m\ 2 m Fz - (m x m z )m Fx + (m y m z )m Fy 
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where V 2 denotes the 2- dimensional Laplacian and a' is the (normalized) damping. The 
dimensionless parameters have roughly the following values: For a permalloy film the 
saturation magnetization is about 640 kA/m. Thus, the Larmor frequency, is about 
2tt ■ 23 GHz, the exchange length, Z ex , is about 6 nanometers, and the dimensionless 
damping, a' is of order 10~ 2 . 

Applying a small-parameter perturbation analysis (rn = am for a a < 1) and 
linearizing Eq. HI one gets 



With the linearized equation, the evolution of excitations at different point contacts can 
be added and patterns studied without integrating the full equations of motion (eq. HJ). 
We present, however, in Sec. [5] computations using the full LLG equations (Eqs. HJ) that 
corroborate this analysis. 

4. Transponder design 

A transponder is a device that can detect spin-wave activity in the propagating medium 
and respond by initiating new spin-waves. It is basically an integrate-and-fire circuit 
[29. 30J, where a storage device (e.g.. a capacitor or an inductor) accumulates charge 
that is discharged rapidly when a certain threshold level is achieved. Such devices are 
widely used in electronics and physics ranging from the macroscopic neon bulb [31] to 
nano-sized memristors [32} [331 EH EHJ [361 EZ] • A basic design for a transponder that we 
have implemented in computer simulations is described here. 

STNOs are used to excite the magnetic moments in thin films and the same (or 
other) contacts can also be used as detectors. The respective alignment of the fixed 
and free layers determines the resistance of the STNO. Thus, small currents (that 
do not excite magnetization dynamics) through the point contacts would serve to 
read the state of the free layer. However, the giant magnetoresistance effect is weak 
when the orientation of one layer tilts only a few degrees, particularly when the layer 
magnetizations are initially collinear. If we include an additional dielectric layer on top 
of the STNO we obtain a magnetic tunnel junction (MTJ), which is either in series 
or parallel with the contact. A ferromagnetic top electrode (with the magnetization 
orthogonal to the free layer) provides the detector with a variable resistance depending 
on the orientation between the magnetic moments of the free layer and the magnetization 
of the ferromagnet on top (see, Fig. [3]^a) and (b)). 

A memristic element together with a capacitor, an inductor and an additional 
resistance are part of the detector as well. Thin films of TiC>2 may be deposited 
on top of the MTJ and may even serve as the dielectric in the MTJ. Memristors 
(memory resistors) [321 [33] are any passive two-terminal circuit elements that maintain 
a functional (nonlinear) relationship between charges and magnetic flux. There are no 
generic memristors. Instead, each device implements a time-varying function of net 




(5) 
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Figure 3. A schematic representation of a tunnel junction on top of the contact of the 
STNO in (a) and next to the excitable point contact in (b). An additional ferromagnet 
layer is required to detect tunnel magnetoresistance (TMR) in both cases. V is the 
applied voltage representing a dc-bias, R is a resistor, V c is the voltage across the 
capacitor, L is an inductor and -RjyjTJ 1S ^ ne resistance of the the tunnel junction. 
In (c), (V — V c )/R provides a current source, and it splits into the current into the 
capacitor (C(dV/dt)) and the current I through the point contact (in the resting 
state the current through the capacitor is 0). In (d), the current source is given by 
(V — V c )/Rmtj- In both cases Rmtj controls the voltage across the memristor and 
whether it stays in the low or high-resistance branch of its I — V characteristics curve. 



charge history. Our circuit will includes additional elements to enable current bias 
bistability in the STNO and the integrate and fire functionality (see Figs. [3]^c) and (d)). 

Point contacts have a very low resistance (R ~ 1 O) and the required current 
densities to excite spin-waves lead to tiny heat dissipation. On the other hand, MTJ 
have resistance of hundreds of ohms and consequently the same current densities would 
produce considerable greater dissipation. Figure Mb) shows a parallel configuration of 
STNO plus MTJ. Additionally, Figs. [3](c) and d depict detecting schemes corresponding 
to both the in-series and parallel configurations of the detectors. 

The process of detection is as follows: 1. A small voltage applied to the detector 
point contact induces a small current (insufficient to produce magnetic excitations in the 
free layer, j < j c ), that is used for detecting the TMR through the tunnel junction. 2. 
The capacitor C is charged and the memristic device is in its high resistance branch. 3. 
When a spin wave perturbation reaches the point contact detector, its resistance changes 
(Rd(0)). 4. The change in resistance produces a voltage change in g with time constant 
(RC). 5. When the threshold in g is reached, there is a jump from the high-resistance 
branch to the low-resistance one and the capacitor, C, is discharged through the point 
contact, thereby producing a pulse excitation. 6. Once the capacitor is discharged the 
resistance at g resets to the high-resistance branch and the charging process restarts. 7. 
During the charging, the value of Rd is no longer important. No detection is possible 




Figure 4. An implementation of a super-threshold memristic circuit in (a). V is 
the applied voltage representing a DC-bias, R is a resistor, V c is the voltage across the 
capacitor, L is an inductor and Rd is the resistance of the point contact plus the tunnel 
junction. (V — V c )/R provides a current source, and it is split into the current into 
the capacitor (C(dV/di)) and the current / through the point contact (in the resting 
state the current through the capacitor is 0). The voltage drop across the memristic 
element is V g = L(dl/dt) + g(I) + IRd- In (b) a simulation of the circuit in (a). 
This figure present the time dependence of both time the current through the point 
contact and the voltage in the capacitor. The resistance Rd is the input signal, and 
it varies in time (bottom panel) and produces discharges of the capacitor through the 
memristic element (it varied three times but only the first two exceeded the threshold 
and produced current discharges). 



until the capacitor is charged. 

The equations describing such detection schemes are 
dV 

C— = (S-V c )/R 1 -I 

(6) 

where R\ and R2 will either be the resistance of the MTJ, Rmtj, or the fixed resistance, 
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R. The function V c = g(I), is the / — ^-characteristic given by the function [31] (see 
Fig. Hb)) 

g(I) = 0.0075/ + 21 exp(-I/2) (7) 

(The I-V characteristics of Pt|Ti02|Pt [3U ESJ [361 E£] or phase change memory cells 
would also work). 

A simulation of the super-threshold memristic circuit is included in Fig. H(c). V 
is the constant applied voltage representing a dc-bias. The voltage at the capacitor, 
V c , together with the current / across the STO are plotted as a function of time. The 
resistance of the MTJ is forced to vary at t — 5, 10 and 15 ns. The forced drops of 
Rmtj cause a shift in the memristor I — V curve and eventually the resistance jumps 
from the high-resistance branch to the lower-resistance one, and at that moment the 
capacitor is discharged through the STNO and the initial state in the memristor (high- 
resistance branch) is rapidly recovered and the capacitor recharged again. Notice that 
small variations of the TMR do not produce a change of state of the memristor and 
consequently the capacitor is not discharged through the STNO (that is the case of the 
variation at t = 15 ms). 

If two transponders (i.e., two arrangements of STNOs) spike together, the excitation 
wave-packets will intersect at different times and positions. These locations and times 
can be determined using the geometry of spin-wave intersections. For example, the 
intersection of two wave fronts will occur on a parabola-like curve, and any transponder 
placed on this curve may be used to detect the event. Notice that delays in the excitation 
of the two spin-wave fronts result in shifts of the spin-waves intersection curves so the 
delays may be detected by differently located transponders. Super-threshold detectors 
would be set to detect only excitations caused by more than one transponder, thus, 
protecting the system against other perturbations —such as thermal noise. Differential 
detection comparing the TMR signal in differently located transponders can also be 
used to determine the spin-wave patterns. 

5. Analysis and Simulation of the Patterning 

We first demonstrate interference patterns between spin waves emanating from two 
sources (two point contacts) 0, [101 E] using numerical simulations of Eq. (@|. Note 
that the waves add as they pass each other, and the intersections trace out parabolic 
curves. Thus, the coincidence of two spin-waves should be detectable by a transponder 
in the ferromagnetic film placed along an intersection parabola. Numerical simulations 
indicate that at the convergence of two waves the magnetization vector is deflected 
approximately twice that when a single wave passes. 

Figure [5] shows spatial dependence of one of the in-plane components of the 
magnetization (\m x \ for convenience), at 3 different times, T\ < r 2 < r 3 , of the waves 
emanating from two STNOs after being excited with current pulses. 

Note that the out-of-plane magnetic moment m z represents the envelope of the 
modulated, in-plane components, m\ = 1 — {m? x + ml) = 1 — |m| 2 , and single contact 
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Figure 5. Spatial dependence of the in-plane component of the magnetization, m x , 
at 3 different times, T\ < T2 < T3. A pair of point contacts are initially excited 
by polarized currents. The excitations diffuse throughout the two-dimensional film. 
Interference patterns are also created in m y . Note that the intersecting waves pass 
through each other, and their amplitudes add. Simulations are performed in a 440 by 
440 nm film, contacts of 20 nm in diameter are separated by 125 nm. The current- 
polarized pulses were switched on for 50 ps and the depicted times correspond to 
T\ = 15 ps, T3 = 175 ps and T3 = 265 ps after the polarized pulse was switched off. 




Figure 6. Three point contacts create magnetic excitations that diffuse and interfere 
throughout the two-dimensional thin film. Intersecting waves create six wave packets, 
and they propagate outwards in specific directions. The upper panels show the 
evolution of the amplitude of one of the in-plane components, \m x \; the lower, show 
the evolution of the perturbation amplitude, |m|. Three different times are plotted, 
i~i < T2 < T3. Simulations are performed in a 440 by 440 nm film, contacts of 20 
nm diameter, and separation between contacts of 65 nm. The current-polarized pulses 
were on for 50 ps and the depicted times correspond to n = 15 ps, t-i = 175 ps and 
T3 = 265 ps after the polarized pulse was switched off. 
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excitations decay radially. In addition to the spatial patterning shown in Fig. that 
diffuses, the frequency, oc |7|-Hq, cause the in-plane components, m x and m y to oscillate. 
Thus, encoding information in the envelope (i.e., magnitude of the oscillation, \m\) 
creates spatial patterns that do not oscillate in time and only diffuse/propagate. Here 
we considered our transponders to be a group of point contacts. 

Figure [6] shows the interference patterns for both the amplitude of one of the 
in-plane components, m x , and the amplitude of the perturbation, \m\, created by 3 
identical current-pulses applied to a triangular array of STNOs. In this particular case 
the excitations interfere in such a way that the spin-wave are enhanced at the vertices 
of an expanding hexagon. The complexity of the in-plane components, m x and m y 
is enormous as they oscillate at a fast rate compared with the diffusion times. The 
envelope, |m|, captures the essence of the perturbation, and indeed reflects those areas 
where excitations are; zeros in the the envelope means no spin-wave activity but zeros 
in either m x or m y carry no information on the intensity of the spin-wave excitation. 



6. Reverberating Structures 

The next step is to design a configuration of transponders in which data can be stored. 
There are several configurations that can maintain stable reverberating activity [1] , and 
hence serve as memory storage units. 
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Figure 7. In (a), differential detection configuration. In (b) and (c) configurations 
requiring converging and intersecting wave-packets. 
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Combinations of transponders (number and location) result in controllable 
interference patterns, but these patterns may be unrecognizable from single-point 
readings. Spatial information may be obtained with detectors in order to distinguish 
whether the detected perturbation was created in, say position a, or position b. We 
note that should single detectors be used, several point contacts will have to fire 
simultaneously in different locations through the film and create constructive patterns 
in the detector that surpass its threshold; which would be set to respond only to higher 
amplitudes larger than that of one excitation. In Fig. [7(a) transponders T\ and T 2 are 
made of two sets of point contacts, one red and one blue. The red one in 7\ spikes as in 
Fig. |6] propagating the perturbation along three directions. The transponder T 2 receives 
a wave-packet in the red detector and nothing in the blue one; only this configuration 
allows T 2 to induce a new excitation that will propagate again through the media and 
be detected by T\. 

Another simple configuration that reverberates is when transponders T% and T 2 in 
Fig. 4b are excited at the same time, then wave-front packets propagate from 7\ and T 2 
intersect at T 3 and T 4 , therefore, the wave-front packets pass through T 3 and T 4 at the 
same time exciting a new pair of waves, which return to T\ and T 2 at the same time, etc. 
The reverberation loop is closed and the activity is periodic with a period depending on 
(D + D')/2, where D is the distance between Ti and T 2 , and D' is the distance between 
T 3 and T 4 . 

Interestingly, one can implement similar reverberating memory using three 
transponders (no differential detection is required here), as in Figure 4c. 

6.1. Simulation of a reverberating structures 

Following is a simulation of a four-transponder reverberating system; four double 
transponders, are arranged in a one dimensional ring as shown in the inset of FigJHl 

Transponders T\ and T 3 initially spike together. The spin-wave perturbations 
diffuse and propagate through transponders T 2 and T 4 where their interference pattern 
will be detected by those transponders (T 2 and T 4 ). The variation of the magnetic 
moments are detected, and after a time delay both transponders (T 2 and T 4 ) spike 
together and create a diffusing spin-wave perturbation that propagates towards T\ and 
T 3 . The cycle is closed once the transponders T\ and T 3 detect the perturbation and 
spike again after a time delay. This reverberation may continue for a long time unless a 
powerful enough noise source appears and interferes in the reverberative configuration. 
Figure M depicts the time evolution of the absolute value of the in-plane magnetization, 
\m\, at transponders 7\ and T 2 for a simulation of eq. HJ 

7. Discussion and Conclusions 

It was shown in Ref. [TJ |3] how reverberating structures may be used to perform 
Boolean and arithmetic computations. Any medium supporting steady progressing 
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Figure 8. Time evolution of the envelope, |m|, at transponders T\ and T2- At 
zero time, T\ and T3 are forced to spike together (positive current during 1 time 
unit) creating a considerable perturbation to the in-plane component of the magnetic 
moment. The perturbation propagates towards the neighbors. The in-plane magnetic 
moment in transponders T2 and T4 increases reaching a maxima after 20 time units. 
Once the incoming wave packet is detected T2 and T4 spike together (positive current 
for 1 time unit) and the pattern keeps repeating. 



waves, vortices, or solitons may be used to implement PWC 

To add more flexibility to the model, we allow for multiple layers of the 
two dimensional structures just described, which may share transponders. Layers 
could communicate via shared transponders and have, at the same time, non-shared 
transponders. In other words, if a transponder on one sheet is excited, it may or may 
not be allowed to produce propagating waves from the equivalent position on a nearby 
layer. Moreover, we may place additional transponders in this second layer wherever we 
choose. One can imagine multiple layers laid on top of each other or arranged in more 
complex topologies. Such designs remain to be investigated. 

An application would be programming transponder arrangements having graded 
reverberating frequencies to create a look-up table, similar to ones in the brain. Consider 
the analogy of a piano keyboard; an array of several devices as in fig [T)^c) having a 
gradient of natural frequencies (i.e., different separations between nodes in reverberators 
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result in different resonant frequencies), are simultaneously exposed to a double-spike 
signal with period r. The device with resonant interspike interval r will respond 
coherently, while the others will not. Namely, the double spike will not match with 
the resonant frequency of others and will create complex reverberating activity from 
which the coherent one can be distinguished. Programming by transponder placement 
is similar to learning in the brain where connections may be strengthen by training. 

Reverberation may be controlled, if desired, by modulating any one of the nodes 
or a common oscillating background signal. It is also possible that a reverberator could 
be interrupted by interference patterns created accidentally from other systems, and 
there may be interference from itself if significant noise or boundary reflections are 
present. For the most part, the computations presented here do not exhibit undesirable 
interference patterns since the calculations are all synchronous and require a threshold 
for recognition, i.e., either two waves converging on a transponder or two transponders 
detecting an interference pattern. Damping of signals, and the introduction of absorbing 
boundary conditions may limit the impact of wavefronts on distant configurations, but 
these aspects as well as the impact of random noise on such reverberators require further 
analysis. 

The PWC methodologies have been shown to enable arithmetic and boolean 
operations and look-up table functionality. Only the latter has been discussed here. 
Other computational aspects of the implementations described in this paper and the 
precision of PWC, which, for example, is limited by the breadth of wavefronts, will be 
the subject of future work. Methodologies shown in this paper focused on spin-waves 
created by STNOs. Other nano-systems may also support wavefronts in service of PWC, 
such are arrays of phase-locked magnetic vortices [38], arrays of Josephson junctions [39], 
and rapid single flux quantum devices [UH ST]. Systems that support solitons of activity 
may also provide useful implementations. 
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Appendix A. A Note on Synchronization and Background Oscillations 

It may be desirable in most applications to provide the system with an external 
background oscillation to ensure the interference occurs between point contacts having 
exactly the same driven frequency and phase. Several modes of oscillation (i.e., several 
distinct frequencies) can be considered [2E]. The background excitations can be created 
either with additional persistent ac-currents applied equally to all contacts or directly 
applying microwave radiation to the whole system of contacts [211 E21 E3J EU [25j EEl EZ] ■ 
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The equation describing the precession amplitude, z(t), of an STNO is reduced to 
the normal form of a Hopf Bifurcation [27J [30] 

^ = (a — \z\ 2 )z + iloz, (A.l) 

where a is now a coefficient involving both positive and negative damping. 

With some mathematical treatment it can be shown there exist intervals around the 
oscillating frequency, u, where the oscillators phase- lock either to an external frequency 
or to another oscillator. 

Appendix A.l. Phase-Locking of STNOs to External Sources 

Consider an STNO forced by an external signal having several frequencies, say < u\ < 
cu 2 < ■ ■ ■ < «m: 

dz M 

— = (a — \z\ 2 )z + i uz + e c n exp(iuj n t). 

dt n =i 
Where e is the combined forcing magnitude relative to other model parameters. The 
complex coefficients in this trigonometric polynomial are 

c n = p„exp(if n ), 

and the variable z may also be written in polar coordinates as 

z = rexp(i9). 

With these notations, the equation for z becomes 

dr M 

— = ( a - r 2 )r + e ^ p n cos(u n t + £ n - 9) 

at 71=1 



and 



— = u + e — sin(w n t + f„ 

at n=l r 



This system has rich phase- locking properties. Let us choose one mode, say ojn and set 

9 = u N t + (j) + ^ N . 

Then 

dr M p n 

— = ( a - r 2 )r + e V — cos((w n - u N ) * + (£„- 6v) - 0) 
dt t^i r 

an< ^ M 

^j- = oj - u N + eJ2~ sin((c; n - u N ) £ + (£„- £at) - </>)• 
at 7i=i r 

Averaging this system over t gives for the second equation 



Spin-wave interference for memory and computation 17 
Therefore, as long as the locking condition 

' < 1 



r 



p N e 

is satisfied where r* is the equilibrium of the radial components, the oscillator will lock 
onto the iV-th mode. 

Appendix A. 2. Mutual Phase-locking of STNOs 

Consider now an array of M STNOs with frequencies, say < uj\ < u 2 < ■ ■ ■ < ojm'- 

dzi M 

—1 = ( a - \zi\ 2 )zi + iUiZi + e V c i>n z n . 

dt n=l 

The variables Zi may be written in polar coordinates as 

Zi = riexp(i9i). 

With this notation, the equations for Zi become 

dn M 

— i = ( a - rj)r-i + e ^ n r n cos(6»i - 6 n ) 

dt n=l 

and 

dB M c r 

— = uJi + e 2^ sin(0 n - 9i). 

at n=1 Ti 

Let us look at one mode, say ojn and assume certain symmetry, = Cj^. Mathematical 
analysis not carried out here shows that the amplitudes of both modes approach the 
same values, r\ — > r*, unless one of the oscillators has zero initial value. In this case the 
difference of phases, A = 9i — 9 N satisfies: 

dA A 

— = Ui - u N - 2ec itN sin A 

As before, we see that the system phase-locks if the condition 



2ec 



LN 



< 1 



is satisfied. 
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